Binder cumulants of an urn model and Ising model above critical dimension 
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Solving numerically master equation for a recently introduced urn model, we show that the 
fourth- and sixth-order cumulants remain constant along an exactly located line of critical points. 
Obtained values are in very good agreement with values predicted by Brezin and Zinn- Justin for the 
Ising model above the critical dimension. At the tricritical point cumulants acquire values which 
also agree with a suitably extended Brezin and Zinn- Justin approach. 



The concept of universality and scale invariance plays a fundamental role in the theory of critical phenomena 
It is well known that at criticality the system is characterized by critical exponents. Calculation of these exponents 
for dimension of the system d lower than the so-called critical dimension d c is a highly nontrivial task 0]. On the 
other hand for d > d c the behaviour of a given system is much simpler and critical exponents take mean-held values 
which are usually simple fractional numbers. 

However, not everything is clearly understood above the critical dimension. One of the examples is the Ising model 
(d c = 4) where despite intensive research serious discrepancies between analytical || and numerical [Q calculations 
still persist. Of particular interest is the value of the Binder cumulant at the critical point. Several years ago Brezin 
and Zinn- Justin (BJ) calculated this quantity using field theory methods f| and only recently numerical simulations 
for the d = 5 model are able to confirm it Q. Some other properties of the Ising model above critical dimension 
are still poorly explained by existing theories. For example, the theoretically predicted leading corrections to the 
susceptibility disagree even up the sign with numerical simulations Q . 

In addition to direct simulations of the nearest-neighbour Ising model, there are also some other ways to study the 
critical point of Ising model above critical dimension. For example, Luijten and Blote used the model with d < 3 but 
with long-range interactions 0. Using such an approach they confirmed with good accuracy the BJ predictions for 
the Binder cumulant. 

In the present paper we propose yet another approach to the problem of cumulants above critical dimension. 
Namely, we calculate fourth- and sixth-order cumulants at the critical point of a recently introduced urn model [^| . 
Albeit structureless, this model exhibits a mean-field Ising-type symmetry breaking. Along an exactly located critical 
line, the obtained values are in a very good agreement with values predicted by BJ. Let us notice that our calculations: 
(i) are not affected by the inaccuracy of the location of the critical point which is a serious problem in the case of 
the Ising model (ii) are based on the numerical solution of the master equation which offers a much better accuracy 
than Monte Carlo simulations. Moreover, we calculate these cumulants at the tricritical point and show that the 
obtained values are also in agreement with suitably extended calculations of BJ. That both the Ising model and the 
(structureless) urn model have the same cumulants is a manifestation of strong universality above the upper critical 
dimension: at the critical point not only the lattice structure but also the lattice itself becomes irrelevant. What really 
matters is the type of symmetry which is broken and since in both cases it is the same Z2 symmetry, the equality of 
cumulants follows. 

Our urn model was motivated by recent experiments on the spatial separation of shaken sand B. In the present 
paper we are not concerned with the relation with granular matter and a more detailed justification of rules of the 
urn model is omitted J8|. The model is defined as follows: N particles are distributed between two urns A and B and 
the number of particles in each urn is denoted as M and N — M, respectively. Particles in a given urn (say A) are 
subject to thermal fluctuations and the temperature T of the urn depends on the number of particles in it as: 

T(x) = T + A(l - x), (1) 

where a; is a fraction of a total number of particles in a given urn and To and A are positive constants. (For urn A 
and B, x — M/N and (N — M)/N, respectively.) Next, we define dynamics of the model ||: 

(i) One of the N particles is selected randomly. 

(ii) With probability exp[^-iy] the selected particle changes urns, where x is the fraction of particles in the urn of a 
selected particle. 

To measure the difference in the occupancy of the urns we define 
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In the steady state the flux of particles changing their positions from A to B equals to the flux from B to A. Since 
the selected particles are uncorrelated, the above requirement can be written as: 

< M > exp[— -1 =< N - M > expf— — ; — — -1, (3) 

VV T{<M/N>y VV T{<(N - M)/N >) J ' KJ 

or equivalently 

(t:+ < e > ) exp [ — - -1 = (i- < e >)exp[— - -1. (4) 

Analysis of eq. (^) shows Q that on the (A,T ) phase diagram symmetric (e = 0) and asymmetric (e ^ 0) solutions 
are separated by the critical line which is given by the following equation 

T = - A/2, 0<A<| (5) 

The critical lines terminates at the tricritical point: A = |,2q = § X ■ Let us notice that a random selection of 
particles implies basically the mean- field nature of this model. Consequently at the critical point [3 — 1/2 and 7 w 1 
(measured from the divergence of the variance of the order parameter), which are ordinary mean- field exponents. 
However, the calculation of the dynamical exponent z gives z = 0.50(1) || while the mean-field value is 2. We do not 
have convincing arguments which would explain such a small value of z. Presumably this fact might be related with 
a structureless nature of our model. 

Defining p(M, t) as the probability that in a given urn (say A) at the time t there are M particles, the evolution of 
the model is described by the following master equation 

p{M, t+l)= N ~^ +1 p(M - 1, t)u(N - M + 1) + — ^p(M + 1, t)u{M + 1) + 
p(M,t){^[l - u(M)] + ^^[1 - «(JV - M)\ } for M = 1,2.. .N - 1 
p(0, t + 1) = ip(l, t)w(l) + p(0, t)[l - w(JV)], 

p(iV, t + 1) = ^p(iV - 1, t)w(l) + p(tf, t) [1 - w(JV)], (6) 

where u>(M) = exp[ T( - A/ y jv ^ ]. Supplementing the above equations with initial conditions one can easily solve them 
numerically 

Cumulants that we calculate are defined as 

^4 = — 20, 2:6 = - n -3 (7) 



where 
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M=0 



<e ">= E(^-2 )np(M,oo) (8) 



and the symbol of infinity indicates that we take the long-time (steady-state) solutions of the master equation (g). 
Calculations are made for A = |, i, i and | and for each A the value of T is calculated from eq. (^|). Thus, the last 
point is the tricritical point and the remaining ones are critical points. Numerical results are presented in Figs. [l]-[|. 

Before discussing our results further, let us briefly describe the BJ approach. To calculate cumulants above the 
critical dimension they used the Ginzburg-Landau- Wilson model. Then, they calculate the effective action restricting 
the expansion only to the homogeneous contributions (the lowest-mode approximation). Since at criticality the 

quadratic (in the order parameter) term vanishes in such an expansion and the leading term is quartic which implies 
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that the probability distribution has the form p(x) ~ e~ x , where x is a rcscalcd order parameter. Calculations of 
moments for such a distribution are then elementary and one obtains 
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XA = i [r( I )]4 ~ 2 - 188440 -' = ~ 6- 565319 -- (9) 

The fact that one can restrict the expansion of the free energy to the lowest order term is by no means obvious || . 
Such a restriction leads to the correct results but only above critical dimension where the model behaves according to 
the mean-field scenario with fluctuations playing negligible role. For d < d c additional terms in the expansion are also 
important and cumulants take different value. Numerical confirmation of the above results requires extensive Monte 
Carlo simulations, and a satisfactory confirmation was obtained only for £4 |p|,^0[. 

Omitting detailed field theory analysis, we can extend the BJ approach to the tricritical point. At such a point also 
the quartic term vanishes which makes the sixth-order term the leading one and the probability distribution gets the 
form p(x) ~ e~ x . Simple calculations for such a distribution yield 
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FIG. 1. The fourth-order cumulant x&(N) as a function of 1/N for (from top) A 
Arrows indicate the BJ results for the critical and the tricritical point. 
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FIG. 2. The same as in Fig. [j] but for the sixth-order cumulant xe(N). 

The BJ results (|9|)-(|i0|) arc indicated by small arrows in Figs. |l|-|[ Even without any extrapolation one can see, 
especially for critical points, a good agreement with our results. Data in Figs, [l]-^] shows strong finite-size corrections. 
To have a better estimations of asymptotic values in the limit N — -> 00 we assume finite size corrections of the form 



X4,e(N) = a; 4 , 6 (oo) + AN' 



(11) 
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The least-square fitting of our finite- iV data to eq. ([fl]) gives x^e (00) which agree with BJ values (||)-([h]) within the 
accuracy better than 0.1%. A better estimation of the correction exponent to is obtained assuming that 2:4,6(00) are 
given by the BJ values. The exponent u> equals then the slope of the date in the logarithmic scale as presented in 
Figs. Q% Our data shows that for the critical(tricritical) point u> — 5(3)- 

Let us notice that leading finite-size corrections to the Binder cumulant in the d = 5 Ising model at the critical point 
are also of the form N~ - 5 (with N being the linear system size) (?]]. Moreover, for the tricritical point but d < d c the 
probability distribution is known to exhibit a three-peak structure which is different than the single-peak form 
p(x) ~ e~ x . 
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FIG. 3. Logarithmic plot of xa(BJ) — X4(N) (+) and x§(BJ) — xq(N) (x) as a function for N for A = 0.5. Dotted straight 
lines have slope 0.5. 
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FIG. 4. Logarithmic plot of X4,{BJ) — xn{N) (+) and xe{BJ) — xq(N) (x) as a function fo TV for A = | (tricritical point). 
Dotted straight lines have slope |. 



In summary, we calculated fourth- and sixth-order cumulants at the critical and tricritical points in an urn model 

which undergoes a symmetry breaking transition. Our results confirm that, as predicted by Brezin and Zinn- Justin, 
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the critical probability distributions of the rescaled order parameter has the form p(x) ~ e~ x . Similarly, for the 
tricritical point our results suggest that p(x) ~ e~ x . 

Although in our opinion convincing, the results are obtained using numerical methods. It would be desirable to 
have analytical arguments for the generation of such probability distributions. It seems that for the presented urn 
model this might be easier than for the Ising-type models. Let us notice that for the simplest urn model, which was 
introduced by Ehrenfest |fl2|| , the steady-state probability distribution can be calculated exactly in the continuum 

limit of the master equation and the result has the form p(x) ~ e~ x , where x is now proportional to the difference 
of occupancy e. In the Ehrenfest model there is no critical point and we expect that a distribution of the type e~ x 
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might characterize our model but off the critical line (in the symmetric phase). We hope that when suitably extended, 
an analytic approach to our model might extract critical and tricritical distributions as well. Such an approach is left 
as a future problem. 
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